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We have developed a track fitting procedure based on Kalman Filter technique for an Iron 
Calorimeter (ICAL) prototype detector when the detector is flushed with single muon tracks. The 
relevant track parameters i.e., momentum, direction and charge are reconstructed and analyzed. 
This paper discusses the design of the prototype detector, its geometry simulation by Geant4, and 
the detector response with the cosmic ray muons. Finally we show the resolution of reconstructed 
momenta and also the charge identification efficiency of and fi~ events in the magnetized ICAL. 



I. INTRODUCTION 



Primary focus of the Iron Calorimeter at India-based 
Neutrino Observatory (INO) [l[ is to study the inter- 
actions involving the atmospheric muon neutrinos and 
anti-neutrinos. The detector will measure the momen- 
tum, direction i.e., tracking, and the sign of the electric 
charge of muons that are produced by the charged cur- 
rent interactions of the muon neutrinos with the detec- 
tor material. As a first step, a prototype of ICAL has 
been already set up at Variable Energy Cyclotron Cen- 
tre (VECC), Kolkata. The prototype will mainly track 
cosmic ray muons. Experience with the prototype will 
be very useful to install much bigger ICAL detector. 

We are involved to develop a track reconstruction tech- 
nique for this prototype detector response. Track recon- 
struction is a process by which one can get an idea of 
the trajectory of the particle inside the detector. This 
track reconstruction process can be split into two meth- 
ods, track finding and track fitting. Different approaches 
to both track finding and the reconstruction of the ini- 
tial track parameters are investigated. In present study 
for the track fitting, "Kalman Filter" [2| technique is uti- 
lized. This method is a very useful one which can be used 
in both track finding and track fitting simultaneously. 




FIG. 1: INO Iron Calorimeter prototype detector at VECC. 



II. PROTOTYPE DETECTOR AND ITS 
SIMULATION 

A. Detector Geometry 

The ICAL prototype detector Q consisting of 13-layers 
of iron plates, each having a dimension of 2.5m x 2.2m, 
and will provide an effective magnetic field of 1.0 m 3 vol- 
ume. Each iron layer is 5.0 cm thick and inter-spaced 
with 5.0 cm gap as shown in Fig. [TJ It is planned that 
a set of 12 Resistive Plate Chambers (RPCs) [4j, having 
dimensions of 1.0m x 1.0m will be placed inside the iron 
layers and these RPCs will be of both glass-type and 
bakelite-type [||. The coils which are passing perpen- 
dicularly through the iron layers will carry approximate 
current of 500 ampere. To reduce the heating effect due 
to this huge amount of current, the coils will be water 
cooled. The calorimeter will be magnetized with a uni- 
form magnetic field of magnitude 1.0 Tesla (approx.) so 
that it can distinguish the /i + and /i - events from their 
opposite bending natures. 

B. Detector Simulation by Geant4 

CERN library based "detector description and simu- 
lation tool" called GEANT 0, is used to develop the 
detector geometry as well as to simulate the prototype 
detector response when the muons pass through it. 

The simulated detector volume consists of iron layers, 
which are stacked along the y-direction and the pick-up 
strips are placed along the x and z-directions inside the 
gap. An uniform magnetic field of magnitude 1 Tesla 
is applied along the x-direction to magnetize the simu- 
lated iron layers. When primary particles i.e., the muons 
are projected along the y-axis, a clear signature of the 
bending is visible in Fig. [21 

C. Results and Discussions 

The prototype detector will not be able to detect at- 
mospheric neutrinos (v) due to its limited size and very 
high cosmic ray muon background on the earth surface. 
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FIG. 2: Simulated prototype detector response while /j, pass- 
ing through it. 
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FIG. 3: Mean energy deposition of muon inside 2 mm thick 
RPC gas. The energy deposition spectrum, simulated by 
GEANT4, has been fitted with Landau spectra. 



So, high energy atmospheric muon flux (in the range 
of 0.5 GeV to 2.0 GeV) on the earth surface is used 
to simulate the detector response. The mean value of 
the energy deposited by fj, in 2.0 mm thick RPC gas is 
evaluated and the value is approximately 1 keV as shown 
in Fig. [3 



past three decades [9(. Now, we would like to discuss 
about the Kalman Filter algorithm and its application. 

A. Kalman Filter Algorithm 

In the Kalman Filter framework, a track is designated 
as a set of parameters, called the Kalman state vector 
(r), which is allowed to change along the particle's path. 
The procedure starts with a certain initial approximation 
r = ro and refines the vector r, consecutively adding more 
measurements. The optimum estimation of the state vec- 
tor is attained after the addition of the last measure- 
ment [ic| . The estimation of r is governed by the linear 
difference equation 

r k+1 = A k r k + i/fe, (1) 

where A k is a linear operator, v k is the process noise 
between (k — l) t h and k t h measurements. 

If the measurement m k linearly depends on r k , then 

m k = H k r k + rj kl (2) 

where H k is the propagator matrix between the measure- 
ment space and Kalman state, and rj k is the measurement 
noise. 

It is assumed that measurement errors rji and the pro- 
cess noise v§ are uncorrelated, unbiased ({r]i) = (fj) = 0), 
and the covariance matrices V k , Q k from these errors are 
expressed as 

<m-r)J >=v u 

■ ' Q.r (3) 

The conventional Kalman Filter (KF) algorithm consists 
of following four stages : : 

1. Initialization Step :: An approximate value of the 
vector ro is chosen to start with. Its covariance matrix 
is set to Co = Inf 2 , where Inf denotes a large positive 
number. 



III. TRACK RECONSTRUCTION FOR THE 
PROTOTYPE DETECTOR 

Pattern Recognition is a technique that studies the op- 
eration and design of systems that recognize patterns in 
data. In data analysis of HEP experiments, this tech- 
nique depends on the type of detector used. In track- 
ing detectors, the signals generated by charged particle 
have to be grouped into track candidates. This is usu- 
ally called track finding or track search. In calorimeters, 
the task of pattern recognition is to group the signals to 
showers and to compute certain properties of the shower. 
The importance and feasibility of track and vertex fit- 
ting [8| in the more complex environment of experiments 
with electronic detectors have been recognized over the 



2. Prediction step :: 

f fe = A k r k -i, 

C k = A k C k ^A T k . (4) 

3. Process noise :: The process noise describes proba- 
bilistic deviations of the state vector r 

r k = fk, 

Ck = C k + Q k . (5) 

4. Filtration step :: In this step r k is updated with the 
new measurement m k to get the optimal estimate of r k 
and its covariance matrix C k 

K k = Cfc-ffJ (v k + H k C k Hl^j , 
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FIG. 4: Comparison of hits from a single muon track before 
and after fit. 
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(6) 



The matrix K k is known as Kalman gain matrix and the 
value x\ is the total \ 2 — deviation from the measure- 
ments mi, ....m k . In brief, in the track fitting algorithm, 
the track parameters are modified and the reconstructed 
parameters are obtained from the optimal state vector. 



B. Implementation of Kalman Filter Algorithm 
and results 



As mentioned earlier, the Kalman Filter algorithm con- 
sists of a succession of alternating prediction and filter 
steps. In this implementation, state vector r = r(x, y, 
dx/dz, dy/dz, q/p) where x and y are for positions, 
dx/dz and dy/dz are their slopes, and q/p is the charge 
to momentum ratio. In the prediction step, the current 
state vector is extrapolated to the next detector plane 
taking into account multiple scattering and energy loss 
of the corresponding particle. In the filtration step, the 
extrapolated state vector is updated by taking a weighted 
mean with the measurement. This means that after each 
prediction step it has to be decided which measurement 
should be used in the subsequent filter step. Convention- 
ally, the measurement which is closest to the prediction 
is selected for inclusion in the filter. The track fitting 
starts from a seed track, for our study we consider only 
single muon tracks and all the tracks are fully-contained 
inside the simulated prototype detector volume. The ini- 
tial state vector components are chosen zero, including 
the initial momentum which reduces biasness. Further 
more, all the diagonal elements of the state vector co- 
variance matrix are assigned to a large value, whereas 
the non-diagonal elements are set to zero. The covari- 
ance matrix takes care of the multiple scatterings and 
energy losses due to the passage of the charged particle 



inside the calorimeter material. At first, the fit proceeds 
towards downstream direction (from first layer to last 
layer) and the parameters are transported to the first 
downstream active plane. This transportation is gov- 
erned by the 4 th order Runge-Kutta method which takes 
care of the bending in magnetic field. A fully-contained 
single particle track is shown in Fig. [5] before and after 
fitting. The track hits will be included in the fit if they 
corresponds to the least % 2 value. This % 2 value is the 
residual between the projected hits and the measured 
hits. During the propagation of the charge particle in- 
side the matter, we have assumed that the particles do 
not decay nor have large scattering angle. We also calcu- 
late process noise covariance matrix during and after each 
propagation. The covariance matrix due to measurement 
noises are set to zero as all the hit points are obtained ex- 
actly from Geant4 simulation. Muon deposited approxi- 
mately 63 MeV energy on an average while traversing one 
5.0cm detector plane perpendicularly, as calculated from 
Bethe-Bloch equation This energy loss is incorpo- 

rated to modify the fifth component of state vector i.e., 
the q/p value as mean energy loss per radiation length. 
When the downstream fit ends, its last updated state 




FIG. 5: Reconstructed momenta for different incident muon 
tracks momenta. 
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FIG. 6: Resolution of reconstructed momenta. 

vector is again used as the starting value for the fit go- 
ing from downstream to upstream. This upstream fit is 
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FIG. 7: Reconstructed momentum distribution for 1 GeV and 
2 GeV muon momenta. 



performed in the same fashion, starting from same initial 
covariance matrix as used during downstream fit. Since, 
Kalman Filter is a recursive algorithm, so we iterate the 
fit process for 5 iterations for better tuning of the fitted 
parameters. 

The linearity of the reconstructed momenta with the 
incident momenta of the reconstructed tracks are shown 
in Fig. Even though the linearity is found to be very 
good, the slope seems to have deviated slightly from 1 
presumably due to partially contained tracks and fluctu- 
ations in energy deposition in the calorimeter material. 
For subsequent analysis, we have worked with calibrated 
reconstructed momenta to make the slope close to unity. 
Performance of the fitting is characterized by the value of 
reconstructed momentum resolution as shown in Fig. [5] 
Now, 



resolution (%) = — x 



100. 



(7) 



Where a is the HWFM of the gaussian and M is the mean 
value obtained from the fit as shown in Fig. [7] For fully- 
contained muon track i.e., the track whose end hit point 
is inside the detector volume, we can estimate its initial 
energy very well. For this prototype, the muon tracks 
are fully contained upto 1.1 GeV, although the program 



is sensitive enough to reconstruct the partially contained 
events as shown in Fig.[5j However, for momenta beyond 
1.2 GeV, the resolutions have been found to get worsened. 
The main reason is that for higher energy, the particles 
go beyond the 13 detector layers. 

The charge identification efficiency is calculated for /u + 
and ji~ events from their opposite bendings in magnetic 
field and it is observed that the efficiency is upto 99 % 
for 2 GeV muon as shown in Fig. [5] The significance of 
this charge separation is that in future this program will 
be able to distinguish neutrino and anti-neutrino events 
which produce yT and [i + respectively by interacting 
with the atoms of detector material. 
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FIG. 8: Charge identification efficiency for various incident 
momenta of muon events. 



IV. SUMMARY 

We have presented a track fitting simulation proce- 
dure for a cosmic ray prototype detector including all the 
physics associated with the passage of a charged particle 
through matter. We have only assumed the single parti- 
cle muon track per event incident on the detector. It is 
found that the program can estimate both fully contained 
and partially contained events with good resolution. We 
believe that by increasing the no of layers, the resolution 
of fitting can be improved. Finally, the procedure can 
prominently distinguish opposite charges of the /i + and 
fj,~ events. 
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